Critical-doping universality for cuprate superconductors: Oxygen 
nuclear-magnetic-resonance investigation of (Ca^Lai_a,)(Bai 75_a.Lao.25+x)Cu30 
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The critical doping levels in cuprates, where the ground state changes its nature (from an an- 
tiferromagnet to a spin glass to superconductor to metal), are not universal. We investigate the 
origin of these critical doping variations by measuring the in-plane oxygen pa hole density in the 
Cu02 layers as a function of the oxygen density y in (Cai,Lai_i:)(Bai.75_a;Lao.25+a;)Cu30y. This 
is done using the oxygen 17 nuclear quadrupole resonance parameter vq. We compare compounds 
with X = 0.1 and 0.4 which have significant critical y variations and find that these variations can 
be explained by a change in the efficiency of hole injection into the p^ orbital. This allows us to 
generate a unified phase diagram for the CLBLCO system across the entire doping range, with no 
adjustable parameters. 



There are several critical doping levels in the phase di- 
agram of the cuprates where the ground state changes 
m. The first critical doping level is when the long 
range antiferromagnetic (AFM) order is destroyed and 
replaced by a spin glass (SG) state; next supercon- 
ductivity (SC) emerges; then the spin glass is de- 
stroyed; and finally, superconductivity is destroyed and 
replaced by a metallic state. These critical levels ex- 
ist in the phase diagram of all cuprates which can be 
doped over a wide range such as La2-2:Sr^Cu04 (LSCO) 
and YBa2Cu30j^ (YBCO), but they vary between com- 
pounds. Several attempts have been made to construct 
a universal phase diagram but thus far only partial dia- 
grams, of only one or two phases, have been achieved 
[2-^]. One particular example is the phase diagram 
of the (Ca::cLai_:r)(Bai.75_a;Lao.25+2;)Cu30j^ (CLBLCO) 
system shown in Fig. [Ija), which includes four different 
families with x = 0.1, 0.2, 0.3 and 0.4 [7]. This phase di- 
agram clearly demonstrates that the critical oxygen den- 
sities y depend on x and thus are not universal. 

The reason for lack of universality is not clear and 
could be one of many. For example, it is possible that the 
doping efficiency of the Cu02 planes is family-dependent. 
A second option is interlayer coupling, which in CLBLCO 
is X dependent Q; it is conceivable that the interlayer 
coupling determines the critical doping. Another possi- 
bility is that two different kinds of holes are formed in 
the Cu02 planes, and only the "mobile holes" partici- 
pate in the SC mechanisrn; perhaps the level of mobility 
varies between families [^, [ifll • Finally, in the t — J model 
the critical doping where the AFM order is destroyed de- 
pends on t/ J [11|; it could be that t/ J varies between 
CLBLCO families. In this work we investigate the ori- 
gin of the critical doping level variation between different 
CLBLCO families by directly measuring the hole density 
in the Cu02 plane of CLBLCO using the oxygen Nuclear 
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FIG. 1. (color online) (Caa:Lai_a:)(Bai.75_j;Lao.25+a;)Cu30j, 
phase diagram (a) showing the superconductivity, Neel, and 
glass critical temperatures Tc, Tn and Tg respectivly as a 
function of x and y Q ■ The arrows mark the oxygen densities 
where Tjv starts to decrease. Fig. (b) is obtained from (a) as 
described in the text. 



strate below. 

Our investigation concentrates on CLBLCO since its 
phase diagram is smooth and systematic. It also does not 
have abrupt features such as the kink in the LSCO SC 
dome at x 



Quadrupole Resonance (NQR) parameter vq. This pa- 
rameter is extracted from Nuclear Magnetic Resonance 
(NMR) experiments and is directly related to the density 
of holes in the in-plane oxygen orbital as we demon- 



0.125 [12|, or a structural phase transition 
and chain ordering as in YBCO. In fact, CLBLCO has a 
tetragonal structure (similar to underdoped YBCO) for 
all values of x and y. Consequently, the layer equivalent 
to the "chain layer" of the YBCO can hold oxygen in 
both the a and b directions. Therefore, for each family. 
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the parameter y varies between 6.4 and 7.25 and con- 
trols the doping level. The crystal quality and the total 
cation charge are also x- independent . These proper- 
ties reduces the number of variable determining T^. Thus, 
understanding superconductivity in this compound could 
shed light on all other cuprates. 

As X increases from x = 0.1 to a; = 0.4 the maximum 
Tc (T™"^) increases from 57 to 81 K; the glass tem- 
perature Tg decreases; and the Neel temperature (Tn) 
of the parent compound increases from 380 to 425 K. 
Another important feature of CLBLCO is that for each 
family there is an oxygen density which marks a tran- 
sition between a constant Tjv and a decreasing Tjv as a 
function of y, as demonstrated in Fig. [Ija). We denote 
this density as yN', for the x = 0.1,0.2,0.3,0.4 families 
yN = 6.69,6.63,6.52,6.43, respectively. 

Due to the systematic behavior of CLBLCO, the four 
different phase diagrams in Fig. [TJa) can be reduced into 
one unified diagram in three steps 0, 1. Extract- 
ing the values of the in-plane AFM coupling J from T/v 
for each family by dividing out the interplane coupling 
contribution. 2. Dividing J, Tg, and Tc of each family 
by T™"^ of that family. 3. Stretching the oxygen den- 
sity axis for each family around its yjv by a factor K{x), 
namely, introducing the quantity K{x){y — yN), where 
K = 0.113, 0.098, 0.079, 0.069, for the x = 0.1 to 0.4 fam- 
ilies respectively. The resulted scaled phase diagram is 
shown in Fig. [lib). The scaling procedure is somewhat 
different from that in previous works where the stretch- 
ing was done around the oxygen density of T™°^ and a 
different set of Ks was used Multiplying all 

the Ks by a numerical factor would yield equally good 
data collapse. Here the Ks are chosen so that optimal 
doping is at K{x){y — y^) = 0.05, for reasons that will 
become clear below. The ratio of K between the x = 0.4 
and a; = 0.1 families is 1 : 1.62. Using the scaling ter- 
minology the question we address experimentally in this 
work is: can K{x) be explained by the in-plane oxygen 
Pa orbital doping efficiency, or is some more exotic ex- 
planation required? 

For our experiments, sintered pellets of CLBLCO with 
different x values were prepared using standard tech- 

Since only 



17[, and then ground into powder. 



niques 

^^O has nuclear spin but its natural abundance is only 
0.038%, the samples were enriched with this isotope. 
In the enrichment process the samples were heated to 
520°C for five days and then cooled to 320°C for five 
more days in enriched oxygen gas with an isotope frac- 
tion of 40 — 50%. In order to obtain different oxygen 
levels some of the enriched samples were later annealed 
in either natural oxygen or nitrogen environments at dif- 
ferent temperatures for 24 hours, and then quenched in 
liquid nitrogen. When possible, Tc was determined with 
a magnetometer as shown in the inset of Fig. [2] We used 
cryogenic SQUID, at a field of 50e in FC conditions. In 
this inset the normalized magnetization of four samples 



with X = 0.1 and y = 7.105,7.055,7.035,7.01 is shown, 
demonstrating the variation in Tc- The oxygen density 
was obtained from the known CLBLCO phase diagram 
(Fig. [T]). For the very underdoped samples, which are 
not superconducting (see samples 1 and 2 in Fig. [T]), the 
oxygen density was determined by iodometric titration 

m. 

In our NMR experiment all the eleven samples empha- 
sized in Fig. [ija) by enlarged symbols were measured. 
We used constant frequency / = 36.525 MHz and tem- 
perature of 110 K (the dependence of vq on temperature 
is within our experimental error [l9|). Measurements 
were performed over a range of external magnetic fields. 
For each field 20480 spin echo sequences were collected. 
The intensity for each field is the integral over the Fourier 
transform of the raw data. 

There are three oxygen sites in CLBLCO: the planar 
O23, apical O4 and "chain" Oi. The NMR line of the 
Oi site is negligible [2^. The apical oxygen does not 
affect the measured NMR lineshape either. We measure 
at 110 K, where its intensity is much smaller than the 
planar oxygen and its line is wider [2^. As a result, 
the spectrum is dominated by only one oxygen nucleus 
instead of three. However, the ^^^La nucleus appears 
in some of the measurements when the amount of ^^O 
in the sample is small. In our working frequency the 
^■^^La has a central transition at 6.1 T, while the ^^O 
central transition is at 6.33 T. Therefore, in samples with 
a low enough concentration of ^^O the lanthanum signal 
changes the NMR line shape at low fields, but not at high 
fields. 

We present in Fig. [2] NMR lines of three samples: (1) 
X = 0.1 close to optimal doping, (2) x = 0.1 underdoped, 
and (3) x — 0.4 very under doped. Their place on the 
phase diagram is indicated in Fig. [Ija). There is a clear 
difference around 6.43 T between lines (1) and (3) but 
lines (2) and (3) are similar. As we explain below, this 
is a consequence of different densities in samples (1) 
and (3), and similar densities in samples (2) and (3). In 
the lower field regime of the underdoped sample (2) the 
lanthanum signal dominates the spectrum. 

A quadrupole nucleus such as ^^O can be viewed as a 
non spherical charge distribution whose energy depends 
on its orientation with respect to the internal electric 
fields. The nuclear Hamiltonian in an external magnetic 
field H is a sum of the usual Zeeman interaction and 
additional quadrupole term, and is given by 



n = -7r^H(l + a)I- 



sQVzz r„ ^2 r2 I f t2 t2\1 

(1) 

where "7 = 5.77 MHz/T is the gyromagnetic ratio, I is 
the nuclear spin operator in the 7 = 5/2 representation, 
a is the shift tensor and eQ is the ^''O quadrupole mo- 
ment. The second term is written in the Electric Field 
Gradient (EEC) system representation where the EEC 
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FIG. 2. (color online) Raw NMR data of three samples 
marked in Fig[TJa): 1 (closed circles) is x=0.1 close to op- 
timal doping (y=7.105), 2 (open circles) is x=0.1 underdoped 
(y=6.9), 3 (closed squares) is x=0.4 underdoped (y=6.79). 
From the high fields data it is clear that the quadrupole fre- 
quency of samples 2 and 3 are almost identical but different 
from sample 1. The dominant contribution in sample 2 at 
lower fields is from the ^^^La nuclei. The inset shows the 
temperature dependence of the normalized magnetic moment 
for the x=0.1 SC samples with y=7.105 (up triangles), 7.055 
(circles), 7.035 (squares) and 7.01 (down triangles). 



tensor is diagonal, Vzz = is the largest (axial) EFG 
eigenvalue, and rj = (V^x — Vyy)/Vzz is the orthorhombic 
EFG asymmetry. 

In the limit of small quadrupole frequency compared 
to ^hH, the energy difference between two nuclear spin 
states AErn^m-i is given by 21| : 



1 



^jhH[l~a,]~huQ[-{3cos'{e)-l) (2) 
-i77sin2(0)cos(2,/))](TO-i) 



where m is the nuclear spin component parallel to the 
external (Zeeman) magnetic field, at = ax + cry + Oz is 
the diagonal term of the chemical and Knight shifts (we 
neglected the ofi-diagonal terms), and Q and are the 
angles between z and the external magnetic field, and 
the quadrupole frequency is defined as 



eQVzz 



4/(2/- 1) 



(3) 



A resonance occurs when the frequency 
/,„ (//, z^Q,CTi,77, 0, 6*) = /\E„i^„i~i/h equals an ap- 
plied frequency /. In powder, all possible orientations 
and line broadening must be taken into account. 
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FIG. 3. (color online) Raw data of an NMR measurement 
(red squares) and a fit of the high field data to Eq. U (blue 
line). At low fields there is a deviation from the raw data due 
to the ^'^^La line. The difference between the raw data and 
the fit is presented in the left inset. The right inset contains 
theoretical plots of samples with two different uq (top lines) 
and two different rj (bottom lines), as explained in the text. 



Therefore, the spectrum is given by 



5/2 °° _M^7 (-Q-o)' 

I{H)^ W{m) dale Uv'^e '^"Q 

m=-3/2 

(4) 



The line broadening in our experiments has typical values 
of Aiyg < 0.2vQ and Aa, = 0.002ai. Finally W{m) 
represents the weights of the different transitions and is 
taken as fit parameter. The effect of the parameter on the 
line shape is demonstrated in the right panel of Fig. [31 
The high field side of the theoretical lines for vq = 0.95, 
1.03 MHz (with ij = 0.33) and of r; = 0.28 ,0.32 (with 
h'Q — 1.02 MHz) are plotted in the right inset. Arrows 
mark the regions in which rj and vq changes have the 
most effect on the spectrum, and enable us to distinguish 
between these parameters. 

An NMR line of the CLBLCO sample with x = 0.4 
and y = 7.1 (solid symbols) is presented in the main 
panel of Fig. |3l The solid line is the best fit of Eq. |4] to 
the data at the high field side. It gives = 0.98 MHz 
and r] = 0.33. These numbers are similar to previous 
measurements of YBCO [23| • The difference between the 
fit and the data on the low field side, which is caused by 
the lanthanum nuclei, is plotted in the left inset. In order 
to obtain VQ{x,y) from all samples we fit Eq. |4]to the 
NMR spectrum of all the measured samples. The results 
are shown in Fig. 01 



4 



The EFG on a planar oxygen site is induced by the 
electrons and nuclei surrounding the oxygen (22| . The 
principal axis of this EFG (z direction) is parallel to the 
copper-oxygen-copper axis [20| . The two main contri- 
butions to Vzz are : I) holes in the oxygen p„ orbital; 
II) holes and nuclei of the atoms surrounding the oxy- 
gen. The first contribution is directly proportional to 
the number of holes created by the doping process. In 
contrast, (II) has a negligible dependence on doping [23| . 
Moreover, the holes in the oxygen are much closer to 
the nucleus and therefore their contribution to the EFG 
is more significant. 

The classical formula of Vzz is given by Vzz = 
/ p(^) '' '^i^ d^i^ where p and r are the charge density 
and distance from the nucleus, respectively. The value of 
Vzz induced by a different charge distribution given by 
p'il^) = ^^p{^~^), where ^ is a constant, is 
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(5) 



The oxygen pa- electronic wave functions of the dif- 
ferent CLBLCO samples differ mostly in their typical 
length-scale. The characteristic length is proportional 
to the unit cell parameter a. Neutron diffraction exper- 
iments show that a changes only by about one percent 
between the different families, and by 0.1 percent within 
a family [ll]. Therefore, using Eq. [5] and Eq. [3] it is 
expected that A{Vzza^) oc A(i/Qa^) c>c Aup^, where rip^ 
is the hole density in the oxygen p^ orbital and A stands 
for the changes induced by the doping process. 

In Fig. |4] we present vqa^ versus oxygen levels y for 
the two families with a; = 0.1 and 0.4 (Samples 1, 2 and 
3 are the same as in Fig. [1] and [2]) . We can clearly see 
that the rate at which vqg^ increases with increasing y 
varies between the two families. The data from the dif- 
ferent families generate two different linear curves. Two 
straight lines are fitted to these datasets with the con- 
straint that the slopes ratio is 1.62, which, as mentioned 
before, is the ratio of ^(0.4) to K{0.1). The measured 
lyga'^ versus y can be explained well by the two lines. 
When we perform a linear fits with no constrains the 
slopes ratio is 2.2±0.65 which is within error bars equal 
to the if 's ratio. This is the main experimental finding 
of this work. 

Since we have demonstrated experimentally that 
K{x) {y — yN) and argued above that 
ip^ , we conclude that 

Anp„ oc K{x) [y - vn) ■ (6) 



A(i/Qa-^) oc 
A(i/Qa3) oc A; 



Therefore, the doping efficiency of the p^r orbital is family 
dependent. We would like to emphasize again that we 
can only quantify the doping efficiency ratio between the 
two different families {Ks ratio), but not their absolute 
value. Hence the proportionality sign in Eq. [6l The set of 
Ks which generate Fig. [IJb) are chosen to give Arip = 
0.05 at optimal doping according to Hasse et al. (23j . 
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FIG. 4. (color online) The translation of the number of oxy- 
gen atoms per unit cell (abscissa) into the number of pa holes 
(ordinate), as extracted from the NMR data. The x=0.1 fam- 
ily is in black circles, and the x=0.4 family is in red squares. 
The ratio between the slopes is equal to the stretching ratio 
between the families in the scaling process shown in Fig. [T] 
(see text). The place on the phase diagram of samples 1 to 
3 is shown in Fig. [ifa), and their raw NMR data is depicted 
in Fig. [21 The inset shows that measurements of the copper 
quadrupole frequencies in CLBLCO give the same slope for 
both families 



Negative values of An^^ represent Cu02 planes which 
are not doped. Another degree of freedom is the number 
of Pa holes for y < y^. Following Ref. [23] again, by 
setting np^{y < y^) = 0.11 we obtain Tip^ — 0.16 at 
optimal doping. 

There were several attempts in the past to find a re- 
lation between the number of holes in the Cu02 plane 
and oxygen level y of CLBLCO. Chmaissem et al. used 
bond valence summation (BVS) calculations based on 
structural parameters determined by neutron diffraction 
24| . Keren et al. measured the in plane ^"^Cu NQR 
parameter ^^vq which is shown in the inset of Fig. |4j 
and Sanna et al. experimented with x-ray fine struc- 
ture (XFS) [l^. BVS has some theoretical arbitrariness 
and is not completely reliable. ^^I'qiy) shows no family 
dependence because ^^i'q is sensitive to charge on the 
apical O4 Per, Cu 3dx2_^2^ 3dz2_r2 and 4s, and 02,3 Pa- 
holes simultaneously [23|; hence it is not an ideal probe 
and a difference in the slopes of ^^i'q{x, y) could not be 
detected within the experimental error bars. Finally, the 
XFS peak is constructed of three contributions [l3| which 
again limit their resolution. Therefore, none of the three 
attempts could find a difference in the doping efficiency 
of the planes within experimental resolution. The oxygen 
NQR has the advantage of measuring directly the depen- 
dence of Per hole density Up^ on the oxygen level y, and, 
indeed, this probe detects variations in doping efficiency. 
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The physical meaning of Fig. |4] and Eq. |6] is that 
the efficiency of the doping process, namely the injec- 
tion of holes into the oxygen p„ orbital, varies between 
the CLBLCO families. Moreover, the scaling procedure 
leading from Fig.[TIa) to Fig. [TJb) can now be fully justi- 
fied: the third step is needed because the oxygen density 
is not the relevant parameter and one must use Up^ in 
the phase diagram. The second step means that T^'^^ is 
determined by the magnetic super-exchange interaction 
energy scale J, namely, J{x)/T^^^{x) is x independent. 
This leads to the unified phase diagram of Fig.[TJb) which 
is obtained with no adjustable parameters. 

This work was funded by the Israeli Science Founda- 
tion. 
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